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NIR/Optical Counterparts of Hotspots in Radio Galaxies 
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ABSTRACT 

We present new high spatial resolution VLT and VLA observations of a sample of nine low- 
power (Pi.4GHz ^ 10^' W/Hz) radio hotspots. Infrared/optical emission is definitely detected 
in four of the nine observed objects, resulting in a detection rate of at least 45%. This emission 
is interpreted as synchrotron radiation from the electrons accelerated in the hot spots. The in- 
tegrated spectra of these hotspots reveal typical break frequencies between 10^ and 10^ GHz, 
two orders of magnitude higher than typically found in high-power hotspots. This supports 
the idea that in low -power hotspots with their relatively low magnetic field strengths electrons 
emit most of their energy at higher frequencies. A simple spectral ageing analysis would im- 
ply that the emitting electrons have been injected into the hotspot volume less than ~ 10^ 
years ago. We discuss possible scenarios to explain the lack of older electrons in the hotspot 
region. In particular, the extended morphology of the NIR/optical emission would suggest that 
efficient re-acceleration mechanisms rejuvenate the electron populations. 

Key words: galaxies: active - jets, radio continuum: galaxies 



1 INTRODUCTION 

Radio hotspots are regions of enhanced radio emission located at 
the end of the radio lobes of powerful radio sources. These remote 
bright regions (usually at several hundred kpc away from the galaxy 
core) mark the "working surface" of supersonic jets in which their 
kinetic energy is dissipated into the acceleration of relativistic par- 
ticles. These relativistic particles are left behind forming the ex- 
tended diffuse radio lobes while the hotspots advance into the am- 
bient medium. As such, hotspots represent the end point of an ener- 
getic jet: how the energy of the central radio source is transported, 
released and maintained in such remote regions are still open ques- 
tions. However, being relatively isolated regions - away from the 
central galaxy core - hotspots are ideal laboratories for testing the 
energetics and the evolution of a radio galaxy. 

The synchrotron spectral energy distribution between the ra- 
dio and optical regime within hotspots, and their radio morphology 
and polarisation suggest that the relativistic electrons are acceler- 
ated at the shocks that form due to the impact of the relativistic jets 
into the ambient medium (e.g. Meisenheimer 2003). In the con- 
text of this scenario the resulting synchrotron spectrum is expected 
to steepen at higher frequencies due to the severe radiative losses 
that electrons experience downstream of the shock region (Heavens 
& Meisenheimer 1987). This is also the reason why optical detec- 
tions of hotspots are still relatively rare although more observations 
have been performed since the first detections at visual wavelengths 
(e.g., Meisenheimer & Roser 1996; Meisenheimer, Yates & Roser 
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1997; Lahteenmaki & Valtaoja 1999; Brunetti et al. 2001; Prieto, 
Brunetti & Mack 2002; Cheung, Kim & Song 2005). 

At higher frequencies, in the X-ray regime, the search for 
counterparts of radio hotspots turned out to be more successful. 
Since the initial X-ray detection of the Cygnus A hotspots by Har- 
ris, Carilli & Perley (1994), many hotspots were found to emit X- 
rays (e.g., Hardcastle 2003). Many of these detections could be ex- 
plained with synchrotron self-Compton (SSC) emission with the 
magnetic field strength close to the equipartition value. A signifi- 
cant number of X-ray emitting hotspots however appears to have 
a larger X-ray intensity than predicted by SSC. This excess could 
be explained by a magnetic field strength much lower than equipar- 
tition, IC emission from the decelerating jet (Georganopoulos & 
Kazanas 2003) or an additional synchrotron component (Hardcas- 
tle et al. 2004). In particular, it has been shown by the latter authors 
that the more luminous hotspots tend to show SSC emission, while 
the weaker radio hotspots emit the extra synchrotron component. 
Both X-ray and optical synchrotron emission from hotspots offer 
the opportunity to study several of their characteristics which have 
a number of fundamental implications for the acceleration of rela- 
tivistic particles in general. They imply either a loss-free transport 
to the hotspot (e.g. Kundt 1996), or an in-situ production of ener- 
getic particles with Lorentz factors > 10^ (e.g. Meisenheimer et al. 
1989). Transport from the nucleus is however limited by the un- 
avoidable inverse Compton losses experienced by the electrons in 
the relativistic jet (Gopal-Krishna et al. 2001). This makes it dif- 
ficult to explain the optical emission from hot spots at sufficiently 
large distance (i.e.^ 150 — 200 kpc) from the nucleus (Brunetti et 
al. 2003) with this scenario. 

In a seminal paper Meisenheimer et al. (1997) discussed the 
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3C105 8439.9 MHz Stokes I 3C105S NIRK-band 




Figure 1. a) 3C 105: entire source, 8.3 GHz, VLA D-array; b) 3C 105S: 8.3 GHz, VLA BC array; c) 3C 105S: superposition of VLA A-array data (contours) 
and ISAAC K-band image (grey), both smoothed to a common resolution of l!'l x l!'l; d) 3C 105S: NIRK-band emission in original resolution (0"43). 



synchrotron spectra of hotspots, from radio to optical, through the 
analysis of a sample of eight hotspots. Most of these hotspots have 
spectra that could in principle be accommodated within the stan- 
dard diffusive shock acceleration model (Bell 1978) on the basis 
of the observed radio photon spectral index of ^0.5 (Sy v^"). 
In the context of this scenario two other characteristics which de- 
scribe the spectral shape once the electrons have been accelerated 
are the break frequency, i.e. the maximum energy of the oldest elec- 
tron population in the hotspot volume, driven by the cooling of the 
electrons in the post-shock region, and the cut-off frequency, gener- 
ated by the competition between acceleration and loss mechanisms 
in the shock region. Observationally, the radiation spectral indices 
typically increase by Aa ~ 0.5 beyond the break frequency, and 
an exponential spectral cut-off develops at the cut-off frequency. 
Brunetti et al. (2003) have shown that the magnetic field strengths 
are the weaker, the higher the break and the cut-off frequencies 
are, with a scaling between break frequency and hotspot magnetic 
field in the form Vt, °= B^^ in agreement with theoretical expecta- 
tions based on the shock-acceleration model. Remarkably, the con- 
sequence is that the weaker the magnetic field strengths are, the 
greater the probability is to find hotspots with counterparts at opti- 
cal wavelengths (Blundell et al. 1999). 

There are however cases of hotspots that do not fit in the 



shock-acceleration scenario, at least in its simplest form. The ex- 
traordinary spectra of some hot spots, lacking any indication of syn- 
chrotron ageing (e.g., 3C390.3, Prieto & Kotilainen 1997; 3C303, 
Meisenheimer et al. 1997), and/or showing optical extension (e.g., 
Pictor A, Thomson, Crane & Mackay 1995; 3C445, Prieto et al. 
2002), call for an efficient re-acceleration mechanism in situ: basi- 
cally, the continuous work of smoothly distributed accelerators at 
different points along the hotspot region. 

This paper is the third in a series of articles reporting the re- 
sults of a project aimed at putting studies of spectral emission from 
hotspots on a broader statistical fundament. In our first paper (Pri- 
eto et al. 2002) we showed the optical images made of the southern 
hotspot of 3C 445. These images revealed bright knots embedded in 
diffuse optical emission distributed along the region of the impact 
of the jet into the intergalactic medium, and suggested that continu- 
ous re-acceleration of electrons (Fermi-II) is at work in this hotspot 
in addition to the prime (Fermi-I) acceleration mechanism. The sec- 
ond paper (Brunetti et al. 2003) focused on the statistical aspects of 
optically detected hotspots identifying the magnetic field strength 
in hotspots as the key parameter which determines the probability 
of their detection at optical wavelengths. The present paper com- 
piles the complete set of optical data from VLT ISAAC and FORS 
observations, complemented by high-resolution radio images on a 
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Figure 2. 3C195, the host galaxy in K-band. The elongated structure in 
east-west direction discussed by Neff et al. (1995) is a superimposed spiral 
galaxy in combination with a stellar-like object. 

sample of candidate objects. In Sect. 2 we describe the selection of 
our sample and introduce the target sources. Sect. 3 gives the de- 
tails of the observations and the data analysis. In Sect. 4 we show 
the results for the individual targets. In Sect. 5 we compare the 
optical magnitudes with the corresponding radio data. A summary 
is presented in Sect. 6. Throughout the paper we use the follow- 
ing cosmological parameters: Hq = 71 kms~'Mpc~^ £1a = 0.73, 
Hm = 0.27. 



2 SELECTION OF THE SAMPLE 

We constructed a sample of hotspots that could be sufficiently 
bright to be detected at the VLT. The source list is selected from 
the complete sample of radio galaxies presented by Tadhunter et al. 
(1993) which contains 54 sources with measured redshifts z < 0.7 
and -25° < 6 < 10° taken from the Wall & Peacock (1985) 2.7- 
GHz sample. All these sources have been visually inspected using 
high-resolution maps published by Morganti, Killeen & Tadhunter 
(1993), Morganti et al. (1999), Kapahi et al. (1998), Leahy et al. 
(1997) and Hardcastle et al. (1998). In view of the expected an- 
ticorrelation of detection probability and magnetic field strength 
(Blundell et al. 1999; Brunetti et al. 2003) we have a priori ex- 
cluded the brightest and most compact radio hotspots which emit 
most of their energy at low frequencies (radio/mm-band). Then we 
adopted a basic synchrotron model for each radio hotspot by as- 
suming a power law spectrum with injection spectral index a — 0.6 
and a radiative break, with Aa = 0.5, at a frequency between the ra- 
dio and the NIR/optical band (e.g. Meisenheimer et al. 1997). We 
finally used the normalization of the synchrotron spectrum given 
by the radio data for each hotspot and assumed a cut-off frequency 
of their synchrotron spectrum at V = 3 • lO'^ Hz. A final sample 
of 10 radio hotspots turned out to be sufficiently bright to favour a 
significant detection in the NIR with VLT. The selected objects are 



up to 3 orders of magnitude less luminous (Tab. 4) than the bright 
and compact hotspots in Cygnus A (PgjGHz ~ lO^^W/Hz). 

For all of them, high resolution (up to 0"2 FWHM) radio maps 
were obtained by us or collected from the literature and are used for 
comparison. In the following we also present maps at lower reso- 
lution to show the large-scale morphological context of the more 
detailed high-resolution maps of the hotspots. We also refer the 
reader to the images presented by Black et al. (1992) [for 3C227 
and 3C403], Leahy et al. (1997) [for 3C 105, 3C327, 3C445] and 
Ncff, Roberts & Hutchings (1995) [for 3C 195] for further infor- 
mation on the overal structure of the sources. 

2.1 The target sample - Radio properties 

2.1.1 3C105 

3C 105 at redshift z = 0.089 is hosted by a Narrow Line Radio 
Galaxy (Tadhunter et al. 1993). The most striking properties of 
this source are the profusions in northern and southern direction of 
the northern and southern lobes, respectively. These could be inter- 
preted as remainders of a previous jet direction. Both hotspots are 
well-defined with the southern one much brighter than its north- 
em counterpart (see also Fig. 7 in Leahy et al. 1997). For our 
near-infrared detection experiment we have selected the southern 
hotspot complex (Fig. 1). Leahy et al. (1997) interpret the two 
northern components of this complex morphology as jet emission. 
In their scenario, the true jet termination hotspot would be the knot 
at a = 04''07'"25;'6 and 5 = 03°40'46", while the southernmost 
knot is supposed to be formed by a recent disconnection event (Cox 
et al. 1991), i.e. through the impact of a spinning jet on ambient ma- 
terial of the cocoon wall, subsequent deceleration and formation of 
a new shock front. The rich substructure of this hotspot is seen as 
an indication of a rapidly evolving, transient system, presumably 
influenced by the superposition of a collimated backflow structure 
(Leahy et al. 1997). 

2.1.2 3C195 

With a redshift of z= 0. 1 1 the Narrow-Line Radio Galaxy (Mor- 
ganti et al. 1997) 3C 195 is the most distant source in our sam- 
ple. Neff et al. (1995) describe it as a large double-lobe source 
with bright hotspots and a complex lobe structure, suggesting small 
changes of the radio axis. The optical protrusion, or jet, orthogonal 
to the radio axis, noted by these authors, turns out to be an underly- 
ing spiral galaxy as can be seen at higher resolution (Fig. 2). Both, 
southern and northern hotspot, were selected as targets for the VLT 
observations (Fig. 3). 

2.1.3 3C227 

3C227 is a Broad-Line Galaxy (Simpson et al. 1996) at redshift 
z = 0.0861. Black et al. (1992) noted that there is no obvious jet 
detection in either lobe. The hotspots in both lobes have a complex 
morphology. The western one consists of two distinct knots (Fig. 
4), the eastern one of those is again composed of two knots which 
remind of a primary - secondary hotspot scenario, similar to the 
situation in 3C 105S. Lonsdale & Barthel (1998) have discussed 
possible scenarios for such double hotspots which can be created 
either via a beam jitter model (as proposed for 3C 105S) or a beam 
deflection model in which the end point of the jet is not marked by 
the primary hotspot but the secondary, which is still being fed by 
the deflected, decelerated beam. On the contrary, 3C 227E shows a 
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Figure 3. a) 3C 195: entire source, 8.5 GHz, VLA B-array; b) 3C 195N: superposition of VLA A-array data (contours) and ISAAC K-band image (grey), 
both smoothed to a common resolution of 0"55 x 0"55; c) 3C 195S: superposition of VLA A-array data (contours) and ISAAC K-band image (grey), both 
smoothed to a common resolution of 0"55 x 0"55. 



more filamentary structure (Fig. 5), with an extended area of max- 
imum emission. Both hotspots were targets of our optical observa- 
tions. 



2.1.4 3C327 

3C327 is a Narrow-Line Radio Galaxy (Simpson et al. 1996) at 
redshift z = 0. 1039 (Fig. 6a). The eastern hotspot complex consists 
of a primary knot and a more diffuse secondary one north of it. 
Both features are embedded in a region of diffuse extended emis- 
sion (Fig. 6b). The dominant feature in the western hotspot is the 
slightly curved structure at the very end of the hotspot complex. 
Leahy et al. (1997) mention the possibility that the brighter struc- 
ture in the south is part of the radio jet (Fig. 6c). We performed 
optical observations of both hotspots. 

2.1.5 3C403 

The Narrow Line Radio Galaxy 3C403 at redshift 0.059 belongs 
to the peculiar class of X-shaped sources which combine two pairs 
of lobes extending in different directions from the same core. One 
of these lobe systems is typically more extended and diffuse, while 
the other one is shorter, brighter and clearly confined. The more 
extended lobes would be the remainders of an earlier phase of jet 



activity, while the shorter, brighter systems would mark the current 
working direction of the jet. Black et al. (1992) show the overall 
morphology (their Fig. 13). Hardcastle et al. (1998) suggest that the 
knot F6 is the primary hotspot in the western lobe. For our detection 
experiment we focused on the outermost western hotspot (Fig. 7) 
at the edge of the younger lobe system which is described as an 
unusually uniform hotspot without apparent substructure by Black 
et al. (1992). 



2.1.6 3C445 

The host galaxy of 3C445 is a Broad-Line Radio Galaxy (Era- 
cleous & Halpern 1994) at a redshift of 0.05623. This source has 
been defined as an example of recurrent activity sources. In addi- 
tion to the bright hotspots in a distance of 315 kpc (north) and 275 
kpc (south) of the core, two additional internal hotspots at 47 kpc 
(north) and 77 kpc (south), well aligned with the outer hotspots and 
the nucleus are found. These are the characteristics of a Double- 
Double Radio Galaxy (Lara et al. 1999, Schoenmakers et al. 2000). 
Both outer hotspots were targets of our study (Figs. 8 & 9). 
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3 OBSERVATIONS AND DATA REDUCTION 

3.1 Near Infrared observations with ISAAC 

The sample of hotspots listed in Table 1 were observed in the near- 
IR with ISAAC at the VLT in the J-, H- and K-bands. An exception 
were the hotspots 3C 227E and W and 3C 327W, which were ob- 
served in K-band only as these are the faintest targets in the sample. 

The pixel scale of ISAAC in those bands is 0"148/pix. All ob- 
servations were done in service mode to garantee optimal weather 
conditions. In each band the achieved spatial resolution was excep- 
tional, with a FWHM < Of'6 on average. Observations were col- 
lected following a random dithering pattern about the science target 
position. Co-addition of the individual frames was done by apply- 
ing the standard shift-and-add technique. The total integration time 
accumulated on the target was on average about 45 minutes per 
band. In Table 1 a compilation of the ISAAC observations is pre- 
sented. New data which were not used in the first two papers are 
marked with an asterisk in Col. 5. 

3.2 Optical observations with FORS 

Follow-up observations of the hotspot sample were carried out with 
FORS at the VLT in the R and B bands, with the exceptions of 
3C 327W and E, 3C 403W, and 3C 445N. 3C 327E is confused by 
a background galaxy, 3C 327W is too faint for optical observations 
and was not detected in the MR; 3C 403W was not detected in the 
NIR either and the follow-up optical observations were done in I- 
band only. Only 3C 445S was in addition observed in the I and U 
bands. The pixel scale provided by the FORS configuration used 
with these bands is Of'20/pix. The observation technique used was 
the same as that described above with ISAAC. The total integration 
time is about 1 h in the I- and B-bands, 0.5 h in the R-band, and 



Table 1. ISAAC observations of the candidate hotspots. Records with an 



asterisk at the observing 


date mark new data, which were not used in 


first two papers. 










source 


band 


FWHM (") 


Magn. 


observing date 


3C 105S 


K 


0.43 


17.8 ±0.2 


Aug 20 2001 


3C 105S 


H 


0.46 


18.3 it 0.1 


Sep 24 2002 


3C 105S 


J 


0.70 


20.2 ±0.2 


Aug 23 2002 


3C 195N 


K 


0.40 


>20.5 


Apr 09 2001 


3C 195N 


H 


0.34 


>21.6 


Jan 21 2003* 


3C 195N 


J 


0.44 


>22.9 


Jan 21 2003* 


3C 195S 


K 


0.40 


21.3±0.2 


Apr 09 2001 


3C195S 


H 


0.34 


>21.6 


Jan 21 2003* 


3C195S 


J 


0.44 


>22.9 


Jan 21 2003* 


3C 227E 


K 


0.47 


20.6 ±0.3 


Apr 18 2002 


3C 227WE 


K 


0.58 


19.8 ±0.3 


Apr 18 2002 


3C327E 


K 


0.37 




May 11 2001 


3C327W 


K 


0.41 


>22.1 


Apr 18 2001 


3C403W 


K 


0.53 


>22.8 


May 11 2003* 


3C403W 


H 


0.56 


>22.3 


Jun06 2001 


3C403W 


J 


0.62 


>23.6 


Jun08 2001 


3C445N 


K 


0.31 


19.9±0.3 


Jun 06 2001 


3C445N 


H 


0.38 


>21.4 


Sep 18 2001 


3C445N 


J 


0.46 


22.1 ±0.3 


Sep 18 2001 


3C445S 


K 


0.75 


19.0 ±0.1 


Jun 07 2001 


3C445S 


H 


0.59 


19.6 ±0.2 


Jul 06 2001 


3C445S 


J 


0.49 


20.2 ±0.1 


Jun 23 2001 



1 .6 h in the U-band. The spatial resolution achieved in these obser- 
vations was on average FWHM ~ 0'.'6. In Table 2 a compilation of 
the FORS observations is presented. New data which were not used 
in the first two papers are marked with an asterisk in Col. 5. 
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3C227 8465 MHz Stokes I 




Figure 5. a) 3C 227: entire source, 8.5 GHz, VLA BC-array; b) 3C 227E: superposition of VLA A-array data (contours) and ISAAC K-band image (grey), 
both smoothed to a common resolution of 0"55 x 0"55; c) 3C 227E: NIR K-band emission in original resolution (0"47). 



o 3C327 8350 MHz Stokes I 




Figure 6. a) 3C 327: entire source, 8.4 GHz, VLA BC-array; b) 3C 327E: superposition of VLA A-array data (contours) and ISAAC K-band image (grey), 
both smoothed to a common resolution of 0"52 x 0"52; c) 3C 327W: superposition of VLA A-array data (contours) and ISAAC K-band image (grey), both 
smoothed to a common resolution of 0"50 x 0"50. 
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Table 2. FORS observations of the candidate hotspots. Records with an 
asterisk at the observing date mark new data, which were not used in the 
first two papers. 



source band FWHM (") Magn. observing date 



3C 105S 


R 


0.6 


22.8 ±0.1 


Nov 26 2003* 


3C 105S 


B 


0.74 


24.5 ±0.3 


Nov 26 2003* 


3C 195N 


R 


0.68 


>25.8 


Dec 18 2003* 


3C 195N 


B 


0.68 


>25.7 


Nov 30 2003* 


3C 195S 


R 


0.68 


25.4±0.2 


Dec 18 2003* 


3C 195S 


B 


0.68 


>25.7 


Nov 30 2003* 


3C 227WE 


R 


0.8 


24.4 ±0.2 


Dec 18 2003* 


3C 227WE 


B 


0.8 


25.4±0.2 


Dec 18 2003* 


3C 227E 


R 


0.8 


23.9±0.2 


Dec 18 2003* 


3C 227E 


B 


0.8 


25.6±0.2 


Dec 18 2003* 


3C403W 


I 


0.65 


>24.0 


Sep 14 2001 


3C445S 


I 


0.66 


21.5±0.1 


Sep 15 2001 


3C445S 


R 


0.58 


22.0±0.2 


Nov 02 2002 


3C445S 


B 


0.56 


22.9 ±0.2 


Nov 02 2002 


3C445S 


U 


0.72 


23.2±0.2 


Jul 03 2003* 



3.3 Radio data 

Between April 2002 and January 2004, we observed all targets, 
with the exception of 3C 227 and 3C 327, with the VLA at 1.4, 5.0 
and/or 8.4 GHz either in A or B configuration, in order to construct 
multi-frequency radio continuum spectra of the hotspots. For the 
source 3C 195 we analysed further archival data at 5.0 and 15.0 
GHz. Information on the observing runs are provided in Table [3] 
The data reduction was carried out following the standard proce- 
dures for the VLA implemented in the NRAO AIPS package. 
For each observing run about 3-5 min, depending on the frequency, 
were spent on the primary calibrator 3C 286. Secondary calibrators, 
chosen to minimize the telescope slewing time were observed for 
1.5 min every 5 min. The flux densities were calibrated according 
to the scale of Baars et al. (1978). At 1.4 GHz it was necessary to 
image a few confusing sources falling into the primary beam in or- 
der to have an accurate measurement of the flux density. The final 
images were produced after a few phase-only self-calibration itera- 
tions. 

In addition, these data were complemented by archive data from 
the NRAO VLA archive, observed as part of the project AB534, 
which were originally published by Black et al. (1992) and Leahy 
et al. (1997). In Table 3 a compilation of the VLA observations is 
presented. New data which were not used in the first two papers are 
marked with an asterisk in Col. 5. 



4 RESULTS 

To identify a possible NIR/optical counterpart of the radio hotspots 
we superimposed the radio images with the NIR/optical data. To 
this purpose, the optical images were brought on the same grid and 
orientation as the radio images. Several sources of the USNO-Bl.O 
survey (Monet et al. 2003) were identified in the images and used 
as reference sources for a transformation of the coordinate systems 
using AIPS task XTRAN. The alignment between radio and optical 
images was cross-checked by superposition of point-like radio and 
optically emitting regions like the host galaxies (when covered by 
the optical image) or other AGN in the field. We estimate the typical 
remaining uncertainties to be of the order of 0"4. In order to assure 
the correct determination of flux densities at the various wavebands 



Table 3. VLA observations of the candidate hotspots 



source 


frequency 


conf. 


beamsize 


beam p. a. 


obs. date 




[GHz] 




[" X "] 


[°] 




3C 105S 


1.4 


A 


1.23x1.18 


35 


Jul 19 2003* 


3C 105S 


4.8 


A 


0.37x0.35 


78 


Jul 19 2003* 


3C 105S 


8.4 


A 


0.54x0.26 


-71 


May 25 1990 


3C 195N, S 


1.4 


B 


6.22x4.01 


-20 


Jan 17 2004* 


3C 195N, S 


4.9 


B 


1.98x1.19 


31 


Aug 18 1986 


3C 195N, S 


8.4 


B 


1.04x0.67 


-15 


Jan 17 2004* 


3C 195N, S 


8.4 


A 


0.32x0.21 


12 


Aug 03 1991 


3C 195N, S 


15.0 


B 


0.50x0.35 


3 


Sep 02 2006* 


3C 227W, E 


8.4 


A 


0.38x0.22 


50 


May 25 1990 


3C 327W 


8.4 


A 


0.35x0.22 


48 


May 25 1990 


3C 403W 


1.4 


A 


1.82x1.37 


58 


Apr 18 2002 


3C403W 


4.8 


A 


0.47x0.35 


55 


Apr 18 2002 


3C403W 


8.4 


B 


0,97x0.85 


-45 


Apr 07 1989 


3C445N 


1.4 


A 


1.43x1.04 


15 


Jul 19 2003* 


3C445N 


4.8 


A 


0.46x0.37 


-23 


Jul 19 2003* 


3C 445N 


8.4 


A 


0.27x0.24 


-31 


Apr 18 2002 


3C445S 


1.4 


A 


1.42x1.02 


11 


Jul 19 2003* 


3C445S 


4.8 


A 


0.46x0.37 


-23 


Jul 19 2003* 


3C445S 


8.4 


A 


0.26x0.22 


12 


Apr 18 2002* 



Table 4. VLA flux densities and angular sizes of the candidate hotspots 



source S1.4 S48 S8.4 S15 Gmaj Omin P8.4GHZ 

[mJy] [mJy] [mjy] [mJy] ["] ["] [10^3 W/Hz] 



3C 105S 


2790 


1268 


873 




7.0 


5.0 


169 


3C 195N 


250 


90 


45 


22 


1.7 


1.4 


13.7 


3C195S 


440 


140 


67 


34 


1.8 


1.1 


20.3 


3C 227E 






47 




1.7 


1.4 


8.5 


3C 227WE 






37 




1.9 


1.2 


3.8 


3C 327W 






21 




2.4 


1.0 


5.6 


3C403W 


94 


50 


33 




5.4 


4.1 


4.7 


3C 445N 


160 


58 


35 




3.5 


2.0 


2.6 


3C445S 


520 


135 


81 




6.0 


3.0 


6.0 



in the same regions, all images were convolved to a common res- 
olution using AIPS task CONVL. The hotspot flux densities and 
angular sizes were measured in the images using the AIPS task JM- 
FIT which performs a Gaussian fit, while in extended structures the 
flux density was obtained by TVSTAT. Observational parameters of 
the hotspots are reported in Tables 1 and 2 for NIR/optical observa- 
tions and Table 4 for the radio observations. The main uncertainties 
of the flux density measurements come from amplitude calibration 
errors, usually between 3% and 5%, while the r.m.s. noise level in 
the image plane is not relevant. 

The optical magnitudes were determined using the IDL-based 
task ATV in various iterations to identify the optimum integration 
area around the targets and correcting for any remaining local back- 
ground instabilities, which however dominate the measurement er- 
rors. Magnitudes were subsequently transformed into jansky using 
standard formulae (Campins, Rieke & Lebofsky 1985; Allen 1972). 
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Figure 7, 3C403: entire source, 8.4 GHz, VLA B-array; b) 3C403W: superposition of VLA B-array data (contours) and ISAAC K-band image (grey), 
smoothed to a common resolution of 1"0 x 1"0. 



3C445 8439.9 MHz Stokes 

a 




" [U \ \ \ L 
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RIGHT ASCENSION (J2000) 
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mjy/beam 
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Figure 8. 3C 445: entire source, 8.4 GHz, VLA C-array ; b) 3C 445N: 8.4 GHz, VLA B-array; c) superposition of VLA A-array data (contours) and ISAAC 
K-band image (grey), both smoothed to a common resolution of 0"35 x 0"35; d) 3C 445N: NIR K-band emission in original resolution (0"31). 



4.1 Optical detections 

We unambiguously detected optical emission connected with the 
radio emission in four of the nine hotspots of the final sample 
(3C327E excluded): 3C105S, 3C227W, 3C445N and 3C445S. 
In addition there are two doubtful cases (3C 195S and 3C 227E), 
for which the successful detection of an NIR counterpart cannot 
be confirmed yet on the basis of the present imaging data. No de- 
tections were found for 3C 195N, 3C327W and 3C403W. The fi- 
nally reached detection rate is 67% (decreasing to 45% if the optical 
counterparts of 3C 195S and 3C 227W are not confirmed). Such a 



success rate is much higher than usually found in the literature and 
may be due to the low magnetic field strengths in these, specifically 
selected, hotspots. 



In Figs. 1 and 3 - 9 we present the NIR K-band emission for 
each of our target superimposed to the radio image, both smoothed 
to the same resolution. In case of a large discrepency between the 
resolutions of the convolved and the original images the optical 
emission is also shown at its original resolution. 
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Figure 9. 3C 445: entire source, 8.4 GHz, VLA C-array; b) 3C 445S: 8.4 GHz, VLA B-array; c) superposition of VLA A-array data (contours) and ISAAC 
K-band image (grey), both smoothed to a common resolution of 0"58 x 0"58. 



4.2 Comments on individual sources 



4.2.1 3C105 



The southern hotspot of 3C 105 is characterized by three well- 
resolved asyinmetric sub-structures, each one with a clear NIR 
counterpart (Fig. 2). The NIR emission at a = 04''07'"25!'22, 
5 = 03°40'51.*'2 does not have a radio counterpart and is considered 
to be an unrelated background source. The northern and southern 
components are located ~ 3"8 (6.2 kpc) and 2"1 (3.4 kpc) respec- 
tively from the central component, which is also the brightest one in 
radio, while in the optical band the southern component accounts 
for the higher flux density. All components have a spectral index 
between 8.4 GHz and K-band of a ~ 1.0. A more detailed analysis 
of this impressive hotspot complex will be presented in a forthcom- 
ing paper. 



4.2.2 3C195 

The high-resolution image of the northern hotspot of 3C 195 re- 
veals the detected NIR emission to be unrelated to the radio emis- 
sion. The situation in the southern hotspot is more difficult to judge. 
The NIR observations show extended emission within the contours 
of the radio emission. However, the shape of the optical emission 
seems to be unrelated to the radio morphology and its barycentre 
is shifted by about one beam size to the southwest with respect to 
the radio peak. This and the relative weakness of the NIR emission 
make a successful detection arguable. 



4.2.3 3C227 

The western lobe is characterized by two asymmetric radio hotspots 
located about 10" (~ 16 kpc) apart from each other. Their flux den- 
sity ratio is Swe^Sww ^ 3.3. Only the eastern component has a 
clear positional counterpart with the NIR emission, while the sit- 
uation of the westernmost component is more uncertain with the 
optical counterpart located outside the compact region of the ra- 
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Figure 10. Examples of model spectral energy distributions of the hotspot 
of 3C 445S for three different break frequencies (assuming ttinj = 0.75) to 
illustrate a comparison between model expectations and data. The models 
are normalised at the radio flux densities. The different break frequencies 
are Vj, = 5.0 • lO"' Hz, Vc/V(, = 1 (solid line); = 1.2 ■ 10'^ Vchh = 50 
(dotted line) and Vfc = 7.0-10" Hz, Vc/vj, = 400 (dashed line). 



dio hotspot (Fig. 4c). The spectral indices between 8.4 GHz and K- 
band are awE ~ 0.9 and aww > 1 .0 for the eastern and the western 
components (assuming that the optical emission is not connected 
with the western radio component), respectively. Hardcastle, Cros- 
ton & Kraft (2007) found X-ray emission from both the eastern and 
western components of the western hotspot of 3C 227. Remarkably 
the radio and X-ray peaks are not coincident and these authors dis- 
cuss possible consequences of the positional offsets on the mecha- 
nisms responsible for the observed emission. Clearly deeper optical 
observations would greatly help to constrain positional offsets be- 
tween the counterparts of the hotspot at different wavelengths and 
to eventually confirm the detection of the optical counterpart of the 
western component of 3C 227W. 

Also the eastern lobe shows underlying NIR emission in the 
area of its hot spot (Fig. 5). These two distinct features are both rel- 
atively weak. The eastern component (Fig. 5c) at the tip of the lobe 
is basically point-like, the elongated feature to the west is fainter 
and roughly coincides with the region of maximum radio bright- 
ness. In principle, this extended NIR emission could derive from 
the radio hotspot, although in order to verify the optical detection 
of the hotspot beyond any doubt much deeper optical imaging is 
required. A comparison with the X-ray emission reported by Hard- 
castle et al. (2007) remarkably shows coincident positions of the 
extended NIR emission with a chain of point-like X-ray emitters 
(Fig. 4 in Hardcastle et al.). 



4.2.4 3C327 

The western hotspot does not show any significant NIR counterpart. 
This results in a lower Umit for the spectral index a > 1.0. The 
image of the eastern hotspot is confused by an underlying spiral 



galaxy. As it is therefore impossible to decide whether the hotspot 
emission reaches up to the NIR band, 3C327E has been removed 
from our sample. 



4.2.5 3C403 

This source at relatively low Galactic latitude is situated in a star- 
rich environment. None of the significant NIR detections within 
the radio contours seems to be related to the radio hotspot. Kraft 
et al. (2005) report extended X-ray emission from the western 
lobe of 3C 403 and suggest inverse Compton scattering of beamed 
IR/optical photons from the active nucleus. On the other hand, these 
authors find several X-ray and optical emitting knots associated 
with the eastern jet, which however are not covered by our data. 



4.2.6 3C445 

Both hotspots of 3C 445 are characterized by extended emission, 
about 2"5 X 1"5 (~2.7 kpc x 1.6 kpc) for the North and about 
7" X 4" (--^7.5 kpc X 4.3 kpc) for the South, detected also in the op- 
tical band. The northern one is elongated in the NE-SW direction 
(Fig. 8c), while the southern hotspot shows an arc-like structure 
(Fig. 9c). Both hotspots have a spectral index between 8.4 GHz and 
K-band of a 0.9. Prieto et al. (2002) underlined the importance 
of 3C 445S that is a clear example of a hotspot whose properties 
are difficult to reproduce with a simple shock-acceleration model. 
The optical images of the hotspot region show an arc-like structure 
with bright compact knots, identified as local particle accelerators, 
embedded into more diffuse optical emission. These knots are in- 
terpreted as the result of a small-scale dentist drill which deposits 
the kinetic energy of the jet in several locations, while the diffuse 
emission may indicate that in-situ particle re-acceleration, due to 
Fermi-II mechanisms, is at work in the whole hotspot region. 

A detailed analysis of these hotspots will be presented in a 
forthcoming paper. 



5 SYNCHROTRON SPECTRA 

In low-power radio hotspots the radiative life-times of the elec- 
tron populations is expected to be longer than in high-power radio 
hotspots, likely due to lower values of the magnetic field (Blun- 
dell et al. 1999; Brunetti et al. 2003). In such a situation, the break 
frequency may occur at IR and optical wavelengths. 

In the context of the shock-acceleration scenario two impor- 
tant frequencies exist in the synchrotron spectrum of hotspots; the 
break frequency, V;,, that marks the maximum frequency of the 
synchrotron emission from the oldest electrons still confined in 
the hotspot volume, and the cut-off frequency, Ve, that marks the 
maximum frequency of the synchrotron radiation emitted by the 
electrons accelerated at the shock (e.g. Meisenheimer et al. 1989; 
Brunetti et al. 2002). Constraining these frequencies provides fun- 
damental implications on the efficiency of the acceleration of par- 
ticles at the shock and on the physical conditions in the shock 
region (e.g. magnetic field strength). Reliable constraints can be 
obtained if the optical data points sample the spectrum at V > V;, 
(e.g., Meisenheimer et al. 1997; Brunetti et al. 2001; Wilson, Young 

6 Shopbell 2001), i.e. in the case of hotspots emitting most of 
their synchrotron radiation at mm-far-IR wavelengths. On the other 
hand, the bulk of the S3mchrotron radiation from the hotspots in our 
sample is emitted at optical wavelengths and thus it is possible to 
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constrain only Vj,'' example is given in Fig. 10, where we show 
thecaseof 3C445S. 

Thus, as Vh/Vc is poorly constrained we decided to fit the 
spectra of the hotspots with a simple single-injection synchrotron 
model (Jaffe & Perola 1973) that is an over-simplification (with 
Vj, = Vc) but still provides a viable approach to obtain statistical 
measures of V;, or upper limits to in the case of optical non- 
detections. More specifically, for 3C 445N and 3C 445S we did not 
consider the flux density at 1 .4 GHz because of the strong contam- 
ination from the lobe. The resulting break frequencies for most of 
the hotspots range between 10^ and 10*' GHz. The upper limit of 
the break frequency obtained for hotspots without optical counter- 
parts is < 5x lO'* GHz. In the case of 3C 327W the availability of 
only two data points did not allow us to perform any synchrotron 
model on their spectrum. Model fits to the spectra of the hotspots 
are shown in Fig. 11, while fit parameters are reported in Table 5. 

From the break frequency v^r [GHz] it is possible to determine 
the radiative age fjaj [Myr] of the electron populations, once the 
magnetic field strength B [/jG] is known, by: 

frad^l610S-3/2— (1) 
\/Vbr(l+z) 

where the magnetic field equivalent to the cosmic microwave 
background has been neglected being much smaller than the typ- 
ical magnetic fields in hotspots. These fields were computed as- 
suming equipartition conditions (Table 5), following the approach 
of Brunetti, Setti & Comastri (1997) assuming a minimum Lorentz 
factor of the relativistic electrons Yn,in = 100, equal proton and elec- 
tron energy densities, and an ellipsoidal geometry with a filling fac- 
tor of unity (i.e. the hotspot volume is fully and homogeneously 
filled by relativistic plasma). 

If we consider the break frequency inferred from the fits we 
find that the radiative ages of the electron populations with detected 
optical emission located in such magnetic fields are of the order of 
10"' years (Table|5}. These relatively short time scales imply either 
that the injection of relativistic plasma from the jet is not continu- 
ous but happens on a period of ~ lO-' years, or that adiabatic ex- 
pansion in the hotspot takes place in ^ 10"' years and the emission 
from older electrons is thus just faded away. 

In the second scenario, however, adiabatic expansion is ex- 
pected to affect the properties of the emitting plasma in a time scale 
that should be significantly larger than a crossing time-scale of the 
hotspot at the speed of sound (or Alfven speed), IhsI^s > IhsI^- 
The crossing time scale should be smaller than (or comparable to) 
the ages measured in our hotspots yielding the following condition: 

3/2 

< 1 (2) 



Table 5. Physical parameters 



B 



eq 



Ihs 
kpc 



that would require either that the magnetic fields in the hotspots are 
substantially smaller than the equipartition field, or that compact 
(sub-kpc), over-pressured, emitting regions, are embedded in the 
hotspot. If these two conditions are not realised the adiabatic losses 
are not fast enough to act within a time-scale of ~ 10"* years and 
the absence of older electrons in the hotspots is not easy to account 
for in the context of a scenario of continuous jet activity. On the 
other hand, the radiative age of the hotspots measured from the syn- 
chrotron break frequency is based on the assumption that electrons, 
accelerated at the shock, cool in the post shock region in absence 
of any additional (efficient) re-acceleration mechanism. This might 
not be the case since the evidence for diffuse synchrotron emission 
on kpc scale in several hotspots in our sample suggests that in situ 
(additional) particle re-acceleration mechanisms may be at work in 



Source 


Vbr 


a 


Beq 


trad 


distance 




[GHz] 




ll^G] 1 


lO' yr] 


[kpc] 


3C 105S 


1.37-10^ 


0.75 


75 


6.4 


278 


3C 195N 


< 2.7-10' 


0.95 


62 


>6.0 


117 


3C 195S 


5.34105 


1.00 


78 


3.0 


127 


3C 227WE 


3.0-105 


0.65 


126 


2.0 


173 


3C 227E 


1.14- 10« 


0.75 


99 


1.5 


169 


3C403W 


<2.95-10'' 


0.55 


38 


>39 


52 


3C 445N 


6.63-105 


0.85 


60 


4.1 


315 


3C 445S 


8.40-105 


0.80 


68 


2.8 


275 



the hotspot volume. In this case the radiative ages estimated from 
the break frequency should represent a lower limit to the real age 
of the electrons. For instance, Prieto et al. (2002) suggested that 
Fermi II mechanisms with a typical electron re-acceleration time- 
scale ~ 10^ years are required to explain the diffuse optical emis- 
sion in the hotspot 3C 445S. If these mechanisms are slightly more 
efficient than previously thought they may contribute to balance the 
radiative losses of the emitting particles, in which case electrons in- 
jected in the hotspot volume more than 10^ years ago may still be 
emitting at optical frequencies. 



6 SUMMARY 

We present new data of a sample of low-power hotspots (Pi ,4ghz ^ 
10^^ W/Hz) observed at several NIR and optical bands, comple- 
mented by high-resolution radio observations. From this investiga- 
tion we found: 

• Of the 9 targeted hotspots at least 4 show a clear detection at 
the NIR K-band, in several cases up to the optical (R, B) bands. An 
additional two hotspots show NIR features which cannot yet safely 
be identified as the optical counterpart of the radio hotspot and will 
need deeper observations. In any case this results in unprecedented 
detection rates between 45% and 67%. 

• This high detection rate is attributed to the high break fre- 
quencies (10^ to 10^ GHz) of the continuum spectra of the de- 
tected sources that emit the bulk of their synchrotron emission at 
optical wavelength. High values of the break frequencies are in- 
deed expected in the class of low-power hotspots (Blundell et al. 
1999; Brunetti et al. 2003). For comparison, typical parameters of 
previously known optical hotspots reported in the literature have 
equipartition fields > lOO^uG and break frequencies in the range of 
10^ to 10"^ GHz. 

• The large values of the synchrotron break frequencies in opti- 
cally detected hotspots in our sample imply that the emitting elec- 
trons have been injected in the hotspot volume less than 10"' years 
ago. A first obvious possibility to explain the lack of older elec- 
trons in the emitting volume is an intermittent activity of the jet 
with a time-scale ~ 10^ years. A second possibility is that the emis- 
sion from older electrons in the hotspot volume is fading away due 
to fast adiabatic losses, in which case the optical emission is ex- 
pected to be produced from compact (sub-kpc) and over-pressured 
knots. Alternatively, in-situ Fermi-II mechanisms, with a particle 
re-acceleration time scale ~ 2 — 5 - 10"' years, partially balance the 
radiative losses of electrons. In this case, the radiative age from the 
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Figure 11. Synchrotron fits to the spectra of the 8 hotspots with sufficient measurements. 3C 327W was excluded as only one radio frequency and one upper 
limit at NIR K-band is known. 



measurement of the synchrotron break frequency should be consid- 
ered a lower limit to the real age of the particles. 

• In almost all cases the detected NIR/optical emission is ex- 
tended between 3 and 8 kpc. It breaks down in multiple substructure 
indicating that we may be over-resolving the emission. The diffuse 
emission suggests that in situ re-acceleration mechanisms may be 
at work in the hotspot volume, as previously pointed out in 3C 445S 
(Meisenheimer et al. 1997, Prieto et al. 2002). This would further 



suggest that the radiative ages measured from the synchrotron break 
differ from the real ages of the emitting particles. 

ACKNOWLEDGEMENTS 

The National Radio Astronomy Observatory is a faciUty of the Na- 
tional Science Fotmdation operated tmder cooperative agreement 



NIR/Optical Observations ofHotspots in Radio Galaxies 13 



Tadhunter C.N., Morganti R., di Serego Alighieri S., Fosbury R.A.E., 

Danziger IJ., 1993, MNRAS, 263, 999 
Thomson R.C., Crane R, Mackay CD., 1995, ApJ, 446, L93 
Wall J. v.. Peacock J.A., 1985, MNRAS, 216, 173 
Wilson A.S., Young A.J., ShopbeU RL., 2001, ApJ, 547, 740 

REFERENCES 

Allen C.W., 1981, 'Astrophysical Quantities', 3rd reprinted edition. The 

Athlone Press, University of London, Great Britain 
Baars J.W.M., Genzel R., Pauliny-Toth U.K., Witzel A., 1977, A&A, 61, 

99 

Bell A.R., 1978, MNRAS, 182, 147 

Black A.R.S., Baum S.A., Leahy J.P, Perley R.A., Riley J.M., Scheuer 
P.A.G., 1992, MNRAS, 256, 186 

Blundell K.M., Rawlings S., Willott C.J., 1999, AJ, 1 17, 677 

Brunetti G., Bondi M., Comastri A., Pedani M., Varano S., Setti G., Hard- 
castle M.J., 2001, ApJ, 561, 157 

Brunetti G., Bondi M., Comastri A., Setti G., 2002, A&A, 381, 795 

Brunetti G., Mack K.-H., Prieto M.A., Varano S., 2003, MNRAS, 345, 40 

Brunetti G., Setti G., Comastri A. 1997, A&A, 325, 898 

Campins H., Rieke G.H., Lebofsky M.J., 1985, AJ, 90, 896 

Cox C.L., Gull S.F., Scheuer PA.G., 1991, MNRAS 252, 558 

Cheung K., Kim C.S., Song J., 2005, ApJ, 628, 104 

Eracleous M., Halpern J.R, 1994, ApJS, 90, 1 

Georganopoulos M., Kazanas D., 2003, ApJ, 594, 27 

Gopal-Krishna, Subramanian P., Wiita P.J., Becker P.A., 2001, A&A, 377, 
827 

Hardcastle M.J., Alexander P, Pooley G.G., RUey J.M., 1998, MNRAS, 
296, 445 

Hardcastle M.J., 2001, A&A, 373, 881 
Hardcastle M.J., 2003, NewAR, 47, 649 

Hardcastle M.J., Harris D.E., Worrall D.M., Birkinshaw M., 2004, ApJ, 
612, 729 

Hardcastle M.J., Croston J.H., Kraft R.R, 2007, ApJ, 669, 893. 
Harris D.E., CariUi C.L., Perley R.A., 1994, Nat, 367, 713 
Heavens A.E, Meisenheimer K., 1987, MNRAS, 225, 335 
Jaffe W.J., Perola G.C., 1973, A&A, 26, 423 

Kapahi V.K., Athreya R.M., van Breugel W., McCarthy P.J., Subrahmanya 

C.R., 1998, ApJS, 118,275 
Kraft R.P, Hardcastle M.J., Worrall D.M., Murray S.S., 2005, ApJ, 622, 

149 

Kundt W., 1996, Lecture Notes in Physics, 471, 1 
Lahteenmaki A., Valtaoja E., 1999, AJ, 117, 1168 

Lara L., Marquez 1., Cotton W.D., Feretti L., Giovannini G., Marcaide J.M., 

Venturi T., 1999, A&A, 348, 699 
Leahy J.P, Black A.R.S., Dennett-Thorpe J., Hardcastle M.M., Komissarov 

S., Perley R.A., Riley J.M., Scheuer PA.G., 1997, MNRAS, 291, 20 
Lonsdale C.J., Barthel PD., 1998, AJ, 115, 895 
Meisenheimer K., Roser H.-J., 1986, Nat, 319, 459 
Meisenheimer K., Roser H.-J., Hiltner PR., Yates M.G., Longair M.S., 

Chini R., Perley R.A., 1989, A&A, 219, 63 
Meisenheimer K., Yates M.G., Roser H.-J., 1997, A&A, 325, 57 
Meisenheimer K., 2003, NewAR, 47, 495 
Monet D.G. et al., 2003, AJ, 125, 984 

Morganti R., Killeen N.E.B., Tadhunter C.N., 1993, MNRAS, 263, 1023 
Morganti R., Oosterloo T.A., Reynolds J.E., Tadhunter C.N., Migenes V., 

1997, MNRAS, 284, 541 
Morganti R., Oosterloo T.A., Tadhunter C.N., Aiudi R., Jones P., ViUar- 

Martin M., 1999, A&AS, 140, 355 
Neff S.G., Roberts L., Hutchings J.B., 1995, ApJS, 99, 349 
Prieto M.A., Kotilainen J.K., 1997, ApJ, 491, L77 
Prieto M.A., Brunetti G., Mack K.-H., 2002, Sci, 298, 193 
Scheuer PA.G., 1995, MNRAS, 277, 331 

Schoenmakers A.P., de Bruyn A.G., Rottgering H.J.A., van der Laan H., 

Kaiser C.R., 2000, MNRAS, 315, 371 
Simpson C, Ward M., Clements D.L., Rawlings S., 1996, MNRAS, 281, 

509 



by As.sociated Universities, Inc. This work has made use of the 
NASA/IPAC Extragalactic Database NED which is operated by the 
JPL, Califomian Institute of Technology, under contract with the 
National Aeronautics and Space Adimnistration. 



